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Abstract 

Background: There are no risk algorithms for the onset of anxiety syndromes at 12 months in primary care. We aimed to 
develop and validate internally a risk algorithm to predict the onset of anxiety syndromes at 12 months. 

Methods: A prospective cohort study with evaluations at baseline, 6 and 12 months. We measured 39 known risk factors 
and used multilevel logistic regression and inverse probability weighting to build the risk algorithm. Our main outcome was 
generalized anxiety, panic and other non-specific anxiety syndromes as measured by the Primary Care Evaluation of IVlental 
Disorders, Patient Health Questionnaire (PRIME-IVID-PHQ). We recruited 3,564 adult primary care attendees without anxiety 
syndromes from 174 family physicians and 32 health centers in 6 Spanish provinces. 

Resu/ts:lhe cumulative 12-month incidence of anxiety syndromes was 12.2%. The predictA-Spain risk algorithm included 
the following predictors of anxiety syndromes: province; sex (female); younger age; taking medicines for anxiety, depression 
or stress; worse physical and mental quality of life {SF-12); dissatisfaction with paid and unpaid work; perception of financial 
strain; and the interactions sex*age, sex*perception of financial strain, and age*dissatisfaction with paid work. The C-index 
was 0.80 (95% confidence interval = 0.78-0.83) and the Hedges' g = 1 .17 (95% confidence interval = 1.04-1.29). The Copas 
shrinkage factor was 0.98 and calibration plots showed an accurate goodness of fit. 

Conclusions: The predictA-Spain risk algorithm is valid to predict anxiety syndromes at 12 months. Although external 
validation is required, the predictA-Spain is available for use as a predictive tool in the prevention of anxiety syndromes in 
primary care. 
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Introduction 

Anxiety disorders occur in 6-12% of the general population [1- 
2], reaching 18.5% in primary care attendees [3]. Of patients with 
an anxiety disorder, 50% have another comorbid psychiatric 
condition such as depression or another anxiety disorder [4] . This 



results in a loss of 805 quality-adjusted life-years annually per 
100,000 primary care patients, thus surpassing chronic physical 
illnesses such as cardiovascular diseases, hypertension, diabetes or 
chronic obstructive pulmonary disease [5]. 

Whilst great progress has been made in the development of 
effective psychological and pharmacological therapies for anxiety 
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Figure 1. Fiow ciiart of patients tKirough tKie predictA-Spain study and numbers becoming anxious. Footnote to Figure 1: DNA: did not 
attend; TO, T6 and T1 2: baseline, 6 and 1 2 months interview. Anxiety syndromes measured by the Primary Care Evaluation of Mental Disorders, Patient 
Health Questionnaire (PRIME-MD-PHQ). ^^-^ 
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[6], these are either not accessed by most persons [7] or they are 
not applied correctly [8]. Accordingly, in addition to improving 
the process of the clinical care of patients with anxiety, 
recommendations suggest developing preventive programs. The 
latter require a clear understanding of the ass{)ciat(;cl risk factors 
[9]. A prospective study of these risk factors, followed over time 
and modeled on large population-based samples, enables the 
development of a risk model for developing an anxiety disorder. 
This can inform individualized preventive therapies based on this 
perception of overall risk [10]. A risk algorithm (PredictA-Europe) 
to predict the onset of generalized anxiety and panic syndromes at 
6 and 24 months sdready exists for European primary care 
attendees [1 1]. However, a prediction period of 6 months may be 
too short if we consider that the symptoms must last at least 6 
months to fulfill the DSM-IV diagnosis of generalized anxiety 
disorder. Furthermore, a period greater than 12 months may 
include people who have developed more than one episode of 
anxiety. We therefore aimed to develop and internally vahdate a 
risk algorithm to predict the onset of anxiety syndromes at 12 
months in primary care attendees. 

Method 

Design and setting 

We undertook a prospective cohort study with evaluations at 
baseline, 6 and 12 months. Although this cohort was originally 
recruited with the aim of developing a risk model for the onset of 
major depression [12], in this analysis we aimed to predict the 
onset of anxiety syndromes. 

The method has been described in detail elsewhere [12]. The 
predictA-Spain study was conducted with the participation of 32 
health centers and 174 family physicians (FPs) distributed 
throughout Spain: Granada in southern Spain; Saragossa and 
La Rioja in northern Spain; Madrid, capital of Spain, situated in 



the center; Las Palmas in the Canary Islands; and Majorca in the 
Balearic Islands. Each health center covers a population of 
15,000-30,000 inhabitants from a geographically defined area. 
The physicians in each health center work as a group, with 
extensive primary care teams. The Spanish National Health 
Service provides free medical cover to 95% of the population. 
Patients can visit their FP as often as they wish without having to 
pay for it, even when they do so for preventive reasons. Each 
patient is assigned to only one FP, who has gatekeeper functions. 
The health centers taking part cover urban and rural settings in 
each province. 

Sampling and exclusion criteria 

In the six Spanish provinces, systematic random samples from 
physician appointment lists were taken at regular intervals of 
between four and six attendees with random starting points for 
each day. The FPs introduced the study to the selected patients 
and requested their permission before contacting the researcher. 
Participants who gave informed consent undertook a research 
interview within two weeks. The study population was recruited 
between October 2005 and February 2006. Exclusion criteria 
were an inability to understand or speak Spanish, severe mental 
disorder (e.g. psychosis, bipolar), dementia or severe neurological/ 
sensor)' illn(;ss, terminal illness, the person was scheduled to be out 
of the city for more than three months during the 1 2 months of 
foUow-up, and persons (representatives) who attended the surgery 
on behalf of the person who had the appointment. 

Variables 

Outcome measure. The outcome of interest was anxiety 
syndromes over the preceding 6 months as defined by the anxiety 

section of the Primary Care Evaluation of Mental Disorders, 
Patient Health Questionnaire (PRIME-MD-PHQ) [13,14]. The 
Spanish version of the PRIME-MD classifies patients who test 
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positive for panic, generalized anxiety and other anxiety disorders 
[14]. The first two diagnoses match the DSM-IV criteria exactly, 
but the third is nonspecific. We used a dichotomous anxiety 
variable to indicate when any of the three diagnoses of anxiety 
were present in a given patient. 

Potential risk factors measure. We selected 39 potential 
risk factors for which there was evidence of reliability and validity 
in the questionnaires used to evaluate them [12]. Baseline 
measurements were made of all the potential risk factors: 

• Socio-demographic factors: (1) age, (2) sex, (3) marital status, 
(4) occupation, (5) employment status, (6) ethnicity, (7) 
nationality, (8) country of birth, (9) educational level, (10) 
income, (11) owner-occupier of their accommodation, (12) 
living alone or with others. 

• Controls, demands and rewards for (13) paid and (14) unpaid 
work, using an adapted version of the job content instrument 
with 7 items each [12,15]. 

• (15) Debt and financial strain by means of three questions with 
Likert responses [16]: 1) General financial strain: "How well 
would you say you are managing financially these days?" (4- 
Likert); 2) Basic financial strain: "How often does it happen 
that you do not have enough money to afford the kind of food 
or clothing you/your family should have?" (5-Likert); and 3) 
Coping with debts: "How much difficulty do you have in 
meeting the payments of household and other biUs?" (6-Likert). 

• (16) Physical and (17) mental well-being, assessed by the 12- 
item Short Form (SF-12) [17,18] and (18) a question on the 
presence of long-standing illness, disability or infirmity. 

• (19) Alcohol misuse, assessed by the Alcohol Use Disorders 
Identification Test (AUDIT) [19-21]. 

• (20) A screen for lifetime depression based on the first two 
questions of the Composite International Diagnostic Interview 
(CIDI) [22]. 

• (21) Lifetime use of recreational drugs (CIDI) [23]. 

• Brief questions on the quality of (22) sexual and (23) emotional 
relationships with a partner, adapted from a standardized 
questionnaire [24]. 

• (24) DSM-IV diagnosis of major depression in the preceding 6 
months using the CIDI [23,25,26]. 

• (25) A question on taking medication for anxiety, depression or 
stress. 

• Childhood experiences of (26) physical, (27) emotional or (28) 
sexual abuse [27]. 

• (29) Nature and strength of spiritual behefs [28]. 

• (30) Presence of serious physical or psychological disorder, or 
substance misuse problems, or any serious disability in persons 
who were close friends or relations of participants. 

• (31) Difficulty getting on with people and maintaining close 
relationships, assessed using questions from a social functioning 
scale [29]. 

• (32) History of serious psychological problems or (33) suicide in 
first-degree relatives [30]. 

• (34) Satisfaction with the neighborhood and (35) perceived 
safety inside/ outside the home using questions from the Health 
Survey for England [3 1] . 

• (36) Threatening events in the preceding 6 months using the 
List of Threatening Experiences Questionnaire [32,33]. 

• (37) Experiences of discrimination in the preceding 6 months 
on grounds of sex, age, ethnicity, appearance, disability, or 
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Table 2. Model to predict drop-out. 







*Odds Ratio 


95% confidence interval 


P 


Province 


Granada 


Saragossa 


1.302 


0.840-2.018 


0.238 


Madrid 


1.868 


1.191-2.930 


0.007 


Logrono 


1.526 


0.988-2.357 


0.057 


Majorca 


3.522 


2.240-5.537 


<0.001 


Las Palmas 


1.711 


0.996-2.936 


0.052 


Sex 


Female 


Male 


1.382 


1.189-1.607 


<0.001 


Age (years) 


0.986 


0.981 -.992 


<0.001 


Education 


Beyond secondary education 


Secondary education 


1.240 


0.960-1 .600 


0.099 


Primary education 


1.328 


1.043-1.692 


0.022 


Incomplete primary education or illiterate 


1.594 


1.178-2.158 


0.002 


Enough money to afford food or clotliing 


Always 


Often 


1.122 


0.912-1.381 


0.275 


Sometimes 


0.982 


0.768-1.256 


0.885 


Seldom 


2.141 


1.015-4.517 


0.046 


Never 


0.650 


0.286-1.479 


0.304 


Discrimination due to age 


No 


Yes 


1.890 


1.169-3.054 


0.009 


Daily smoker 


Non or ex-smoker 


Up to 10 cigarettes 


1.199 


0.941-1.528 


0.142 


From 11 to 20 cigarettes 


1.316 


1 .034-1 .676 


0.026 


From 21 to 30 cigarettes 


0.873 


0.573-1.331 


0.529 


More than 30 cigarettes 


1.391 


0.773-2.502 


0.271 


Lifetime depression 


No 


Yes 


1.22 


1.048-1.415 


0.010 



(*) Multi-level logistic regression with health center as a random component. We selected variables included in the final model from the 39 measured in this study using 
a threshold for inclusion of p<0.20 in bivariate regression. From the model thus obtained, those variables with p>0.05 were extracted step by step to obtain a more 
parsimonious model. 
doi:10.1371/journal.pone.0106370.t002 



sexual orientation, using questions from a European study 
[34]. 

• (38) Adequacy of social support from family and friends [35]. 

• (39) Two questions about smoking habits [36] . 

Statistical Analysis 

Participants with missing anxiety diagnoses at both follow-up 
points (at 6 and 1 2 months) were excluded. We also excluded those 
with missing anxiety diagnoses at one foUow-up point who had no 
anxiety syndromes at the other. However, we included patients 
with anxiety syndromes at one foUow-up point and missing 
diagnoses at the other (at 6 or 1 2 months), as they met the outcome 



criterion of anxiety at some point over the 12 months. We 
performed multilevel logistic regression with cumulative anxiety 
incidence as the dependent variable and health center as a random 
component. To test the hierarchical data structure we used the 
likeliliood-ratio test of the null model taking cumulative incidence 
of anxiety syndromes at 1 2 months as the dependent variable and 
health center as a random factor versus usual logistic regression 
[Chi2 = 28.94; p<0.0001]. The Intraclass Correlation CoeflTicient 
for Health Center was 0.082 (95% Confidence Interval: 0.039- 
0.166). The likelihood-ratio test of the null model with the variable 
doctor as a random factor versus usual logistic regression was also 
significant [Chi2 = 12.81; p = 0.0002]. The Intraclass Correlation 
Coefficient of the variable doctor was 0.091 (95% Confidence 
Interval: 0.044-0. 180). We then checked the likelihood-ratio test of 
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O E 



the null model with health center and doctor as random factors 
versus the nuU model with only health center [Chi2 = 0.71; 
p = 0.2002]. We therefore decided to use multilevel logistic 
regression with health center as the random component. 

We selected variables using a threshold for inclusion of p<0.20 
to ensure that information lost as a result of exclusion of a variable 
from the equation was minimal [37]. From the model thus 
obtained, those variables with p>0.05 were extracted step by step 
to obtain a more parsimonious model. The usefulness of including 
first-degree interactions was considered, especially the interaction 
age*sex because it has been found previously in depression [38,39] 
and anxiety [40] , as well as the dififerent combinations of age and 
sex with the other variables included in the model. We decided to 
include an interaction in the model when the likelihood ratio test 
was significant at p<0.05. We used inverse probability weighting 
[4 1 ,42] to adjust for a possible attrition bias due to participants lost 
to follow-up. AH reported P values were two-sided. 

The ability to distinguish those who would develop anxiety 
° syndromes from those who would not was assessed using the C- 

S index [43]. To compare the discriminative vahdity between two 

risk algorithms we performed the test for two correlated C-index. 
S' Prediction models derived with multivariable regression analysis 

J are known for overestimating regression coefficients. We used a 

calculation proposed by Copas [44] to estimate overfitting of our 
g. prediction models. We calculated effect sizes using Hedges' g [45]. 

V! Calibration, which is the agreement between the observed 

proportions of major depressive disorder and the predicted risks, 
was studied with calibration plots taking deciles of risk. We 
conducted all analyses using STATA, release 12 [46]. 



E "3 



S % Ethics Statement 

o gj. The predictA-Spain study has been conducted according to the 

# principles expressed in the Declaration of Helsinki. This study 

2 S complies with the Code of Ethics of the World Medical 
Association and was approved by the ethics committees: Ethics 
Committee on Human Research of the University of Granada, 

g 7 J Ethics and Research Committee of Primary Health District of 
Malaga, Ethics Committee for Clinical Research of Aragon 

Q ^ T; (CEICA). Research assistants explained to patients the predictA- 

Q. X Spain study in detail, their commitments and rights, and answered 

g oj ai all questions that patients wanted to ask. All the participants read 

=5 ^ ^ an information sheet and signed consent forms to take part in the 

J! !t (u Study. 

lu .C:! 
■£ c " 

g^-l Results 
1 " ^ 

E g Of the 6,299 primary care attendees approached, 1,251 (19.9%) 

f! S3 . were excluded: 506 (8.03%) were outside the age range (18-75 

d 2 § years); 446 (7.1 %) were either representatives of patients or did not 

^ ^ d -D attend the appointment; 156 (2.5%) had severe mental disorder, 

"i S a? "I dementia or severe neurological/sensory iUness; 63 (1.0%) 

0) si 2 terminal illness; 47 (0.75%) trouble communicating in Spanish; 

5 ^ — ilj m and 33 (0.52%) were scheduled to be out of the city for longer than 

i d d <!' d three months during the 12 months of foUow-up. Of the remaining 

o S J" 3 5,048 patients asked to take part in the study 4,166 (82.5%) gave 

■S 1 ? '5 2 their consent. These were then interviewed at baseline, but 585 

.£ VJ^*:^ ^ ^ (14.0%) had a positive diagnosis of anxiety syndrome (by PRIME- 

iii^i =5 MD-PHQ) and 17 (0.41%) had a missing diagnosis, so they were 

5 > > 5 E also excluded. Thus, our at-risk population comprised 3,564 

.1 .1 ra ^ patients (Figure 1) 



(TJ 

-§ -S .E S ^ Patients' socio-demographic characteristics are shown in 

^llsEd Table 1. 

sJqu| Of the 3,564 patients, 2,420 (68%) were interviewed at 6 

months and 2,294 (64.4%) at 12 months. The variables associated 
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Table 4. Comparison of the base model with the models 
including each of the interactions tested. 



Interactions 


Likelihood-ratios between model without 
interactions and model with each interaction 




chi2 


P 


Sex*province^ 


Sex*age 


7.12 


0.008 


Sex*financial strain 


12.98 


0.005 


Sex*taking 
psychotropic drugs 


2.60 


0.107 


Sex*physical health 


0.00 


0.945 


Sex*mental health 


1.31 


0.253 


Sex*paid work 


2.70 


0.441 


Sex*unpaid work 


1.42 


0.702 


Age*province 


4.24 


0.515 


Age*financial strain 


3.14 


0.371 


Age*taklng 
psychotropic 


0.02 


0.889 


Age'physlcal health 


0.23 


0.630 


Age*mental health 


2.12 


0.145 


Age*paid work 


11.01 


0.012 


Age*unpaid work 


1.36 


0.716 



^he model sex*province did not converge. 
doi:l 0.1 371 /journal.pone.CI 06370.t004 



with drop-outs were province (Madrid and Majorca), sex (male), 
lower age, lower educational level, seldom having enough money 
to afford food or clothing (basic financial strain), discrimination 
due to age, daily smoker, and lifetime depression (see Table 2). 
There were only nine missing values for predictor variables in the 
final models to predict the onset of anxiety syndromes. The 
cumulative 12-month incidence of anxiety syndromes was 12.2% 
(panic disorder 4.5%, generalized anxiety disorder 4.1%, other 
anxiety disorders 5%). 

The final model to predict anxiety syndromes included 9 
variables and 3 interactions (Table 3): province; sex (female); lower 
age; taking medicines for anxiety, depression or stress; worse 
physical and mental life quality (SF-12); dissatisfaction with paid 
and unpaid work; perceived general financial strain; and the 
interactions sex*age, sex*perceived financial strain, and age*dissa- 
tisfaction with paid work. From the 15 interactions tested we 
selected the three that had a p<0.05. Firstly, we compared the 
nine-variable model plus the interaction sex*age and the nine- 
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Figure 3. Incidence of anxiety by sex and level of financial 
strain. 

doi:10.1 371/journal.pone.01 06370.g003 

variable model plus the interactions sex*age and sex*financial 
strain \chi2 = 13.90; p = 0.003]. Secondly, we compared the nine- 
variable model plus the interaction sex*age with the nine-variable 
model plus the interactions sex*age and age*paid work 
[chi2 — 1 1.55; p — 0.009]. Thirdly, we compared the nine-variable 
model plus the interactions sex*age and sex*financial strain with 
the nine-variable model plus the interactions sex*age, sex*financial 
problems, and age*paid work [chi2 = \ I J 6; ^ = 0.008]. The 
likelihood ratios for the interactions analyzed are described in 
Table 4 and charts of the three interactions are shown in 
Figures 2-3-4. 

Concerning the interaction sex*age, whereas the tendency for 
the incidence of anxiety in women could be drawn with a more or 
less descending line with effect from the age of 25 years, the 
tendency for the men rose from the age of 25 years up to 54 years, 
with the peak between 45-54 years (Figure 2). 

The interaction sex*financial strain showed that in men with 
many financial problems the risk for the incidence of anxiety was 
increased four-fold (Figure 3). 

The interaction age*dissatisfaction at work had an antagonistic 
effoct, with dissatisfaction at work in the younger patients 
increasing the incidence of anxiety (Figure 4). 

The shrinkage factor was 0.9793 (shrinkage = 1 indicates that 
there is no overestimation). The C-index and effect size (Hedges' g) 
were 0.80 [95% confidence interval (CI) = 0.78-0.83] and 1.17 
(95%CI = 1.04— 1.29), respectively, and the C-index data are 
represented in Figure 5. There was a slight improvement in the 
C-index in the model with interactions (0.8024; 95%CI = 0.7763- 
0.8284) versus the model without interactions (0.7917; 
95%CI = 0.7653-0.8181); Chi^ = 5.1 1(DF = 1) and p = 0.0238 
(Figure 5). 
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Figure 2. Incidence of anxiety by age and sex. 

doi:10.1371/journal.pone.0106370.g002 



Figure 4. Incidence of anxiety by age and dissatisfaction at 
work. 

doi:1 0.1 371/journal.pone.01 06370.g004 
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Figure 5. Area under the ROC curve for the models with/ 
without interactions. Footnote toFigure 5: predictA: model without 
interactions predictAS: model with 3 interactions. 
doi:1 0.1 371/journal.pone.01 06370.g005 

We incorporated depression (measured by tlie CIDI at baseline) 
as a risk factor candidate to be included in the model to predict the 
onset of anxiety syndromes but, although it was statistically 
significant in bivariate analysis (OR = 2.83; 95% CI = 1.70-4.71; 
p<0.0001), it did not remain in the final model after adjustment. If 
we removed from the final model the variables "worse mental life 
quality (SF- 1 2-mental) and taking medication for anxiety, depres- 
sion or stress" and then included major depression (CIDI), the 
latter reached statistical significance (OR= 1.88; 95%CI= 1.18- 
2.99; p = 0.008); but if we just removed only one of them (SF12- 
mental or taking psychotropic agents) then major depression was 
not significant. Our final model was clearly more discriminant (C- 
index = 0.80; 95%CI = 0.78-0.83) than the model excluding SF12- 
mental and taking psychotropic agents and including major 
depression (C-index = 0.77; 95%CI = 0.74-0.79); with the test for 
its difference being statistically significant [chi2(l) = 14.54; 
P = 0.0001]. 

The calibration showed an accurate goodness of fit (Figure 6). 
The predicted probability cut-point of 10% was associated with 
the greatest Youdeii's J statistic J= Sensitivity + Specificity -1), 
which had good sensitivity (80.2%) but poor specificity (66.5%). 
The predicted probability of 1 1 % reached a specificity of 70% and 
a sensitivity of 75%; while the predicted probability of 13% had 
more specificity (75%) than sensitivity (approximately 69%), see 
Table 5. 

Differences (variations in coefficients >10%) were found 
between the final model and the same model weighted for the 
inverse probability of remaining in the follow-up to 12 months 
(Table 6). Examples of the kinds of participants scoring at 
increasing levels of predicted probability of anxiety syndromes 
are shown in Table 7. 

Discussion 

The predictA-Spain risk algorithm is valid. To our knowledge, it 
is the first algorithm that has been developed and internally 
vahdated to predict the onset of anxiety syndromes at 1 2 months 
in primary care attendees. From the shrinkage factor obtained, it 
can be deduced that coefficients were minimally over-estimated in 
the internal validation process. However, external validations are 
required to apply this risk algorithm in different populations. We 
used multi-level regression because of the hierarchical structure of 
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Figure 6. Calibration plots (mean predicted probability versus 
observed probability of anxiety within deciles of predicted 
risk) of the predictA-Spain risk algorithm. 

doi:10.1 371/journal.pone.01 06370.g006 

the data. This approach improves the accuracy of estimates of 
coefficients and standard errors [47]. Our large sample size and 
the number of events (people developing anxiety syndromes) per 
variable included in the model reduced the risk of selecting 
unimportant variables and failing to include important variables 
[48]. However, the sample size was not large enough to address 
external validation in this study, derivating the algorithm in some 
provinces and validating in the others. Though our sample 
possibly under-represented patients who attend infrequendy [49], 
frequent attendees are more likely to suffer from anxiety disorders 
[50] and therefore are most in need of prevention. 

Important differences were seen between the predictA-Spain 
models with and without inverse probabihty weighting, indicating 
that loss to follow-up might lead to attrition bias and that this 
strategy could provide unbiased estimates of coefficients, even in 
the presence of attrition bias [41]. 

The questionnaire used to evaluate our outcome, PRIME-MD, 
has good rehabUity and vahdity indices [14], but we cannot rule 
out classification bias. Moreover, we only considered generalized 
anxiety and panic disorders as defined by the PRIME-MD. Our 
data do not concern other anxiety disorders such as post-traumatic 
stress, obsessive compulsive or phobic anxiety disorders. 

In terms of C-index and effect size (Hedges' g), the predictA- 
Spain risk score compares favorably with the predictA-Europe, 
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Table 5. Predicted probability cut-points to predict anxiety syndromes at 12 months and tlieir associated validity. 



Predicted probability cut-points 


Sensitivity (%) 


Specificity (%) 


LR+ 


LR- 


>0.0817 


84.8 


59.9 


2.12 


0.25 


>0.0928 


82.5 


64.1 


2.30 


0.27 


>0.0998* 


80.2 


66.5 


2.39 


0.30 


>0.1053 


77.4 


67.8 


2.41 


0.33 


>0.1108 


74.7 


69.8 


2.48 


0.36 


>0.1147 


72.4 


71.0 


2.50 


0.40 


>0.1204 


70.0 


72.9 


2.60 


0.41 


>0.1291 


68.5 


75.1 


2.75 


0.42 


2:0.1442 


65.4 


80.0 


2.96 


0.44 


>0.1699 


60.3 


82.7 


3.48 


0.48 


>0.1860 


55.3 


84.7 


3.62 


0.53 



LR+ Likelihood ratio of the positive test; LR- Likelihood ratio of the negative test. 
*Cutpoint where Youden's J statistic (J = Sensitivity + Specificity -1} was greater. 
doi:1 0.1 371 /journal.pone.Ol 06370.t005 



although the difTerences between the two risk scores might be 
explained in part by the greater homogeneity of the sample in 
Spain. The predictA-Spain is also similar to other risk algorithms 
for the onset of major depression [39,51] as well as for risk indices 
for cardiovascular events [52]. 

The predictA-Spain risk algorithm shared some risk factors with 
the predictA-Europe [1 1] (sex, age, physical and mental quality of 
life, and dissatisfaction with paid and unpaid work) but differed in 
others: (1) taking medicines for anxiety, depression or stress, (2) 
perceived financial problems, and (3) the interactions. Moreover, 
the variable dissatisfaction with paid and unpaid work was 
measured differentiy (a dichotomous combined variable in the 
predictA-Europe and two separate scales in the predictA-Spain). 
The first difference (1) might, in theory, be because that variable 
was not included as a candidate in the stepwise selection process in 
predictA-Europe. 

We are quite sure that our population at risk had no anxiety 
syndromes at baseline, at least during the preceding 6 months as 
defined by the anxiety section of the PRIME-MD. However, we 
cannot rule out that some patients had suffered anxiety syndromes 
prior to that date. In addition, our population at risk could have 
had depression or dysthymia at baseline or before. For any of these 
reasons, patients might be taking (appropriately or inappropriate- 
ly) antidepressants and/ or anxiolytics at baseline. There may also 
be patients (without diagnostic criteria for anxiety or depression 
disorders) who had insomnia and for whom a doctor could have 
prescribed pills to sleep. In the case of taking anxiolytics over a 
long time, along with insomnia, a percentage of the patients could 
also have an addiction disorder. Moreover, the question is phrased 
in such a way that it might include those taking anxiolytics- 
antidepressants or other medicines (vitamins, placebos, etc.), often 
inadequately, for minor emotional problems. 

Including this variable in the final model does not mean that 
taking psychotropic drugs causes anxiety syndromes. If our 
objective is to obtain a risk algorithm, all variables potentially 
associated with the outcome, not necessarily causally, can be 
considered in order to predict as accurately as possible [53]. We 
believe that when patients respond that they are taking medication 
for anxiety, depression or stress, it should be understood as an 
intermediate variable related to one or more mental health 
disorder, such as depression, dysthymia, insomnia, addictions, 
personality disorders, adjustment disorders or other minor 



emotional problems. We have checked this with our data and 
this may be true at least in the case of major depression (see the 
results section); additionally, suffering other mental disorders is a 
well known risk factor for anxiety disorders [54]. 

One hypothesis to explain the inclusion of this variable in the 
predictA-Sparn risk algorithm might be that Spanish patients have 
a tendency to ask their FPs for more psychotropic drugs for 
emotional problems encountered in everyday life and Spanish FPs 
tend to give them more medication. There is some support for this 
hypothesis because Spain is among those European countries 
having a higher use of psychotropic drugs [55]. 

The question "taking medication for anxiety, depression or 
stress" is a litde ambiguous, but it is a good and independent 
predictor of the onset of anxiety syndromes and it is also very easy 
to obtain an answer. 

The second difference (2) may be related with the fact that the 
economic indicators for Spain are below those of the European 
mean [56]. Finally, (3) the fact that the interactions were not 
included in the predictA-Europe might be explained because only 
the interactions between time and variables were explored in the 
model [11]. 

The selection of risk factors was performed with the arm of 
building a predictive model for major depression. However, it is 
known that anxiety and depression share most of their risk factors 
[11,39,51,57]. The predictA-Spain shared seven risk factors with 
the predictD-Spain [39] risk algorithm to predict the onset of 
major depression in primary care, although their coefficients were 
different. However, five risk factors were only included in the 
predictD-Spain, whilst a further four were included in the 
predictA-Sparn. We highlight the inclusion of dissatisfaction with 
paid work and the perception of financial strain, which were good 
predictors of anxiety but not depression [39]. Furthermore, the 
perception of many financial difficulties by men (interaction 
financial strain*sex) quadrupled the incidence of anxiety at 12 
months, reaching 80% (see Figure 3); which might trigger a 
priority prevention for these cases. This may possibly be explained 
by the predominant role of the male in Spanish culture to 
maintain the family financially. Thus, given a situation of 
maximum hardship the male could either be affected more or 
have less capacity than the female to cope adequately. Concerning 
the antagonistic interaction sex*age (see Figure 2), a similar result 
was found in a study of the prevalence of anxiety in Europe [58] . 
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Table 6. Weighted and unweighted models* 


by the Inverse 


probability of remaining 


In the follow-up 


to 12 months. 






Model weighted 




Model unweighted 




Risk factors 


B 


S.E. 


P 


a 


S.E. 


P 


Constant Province 


0.475 


1.578 


0.161 


0265 


0.893 


0.767 


Granada 


Saragossa 


1.064 


0.241 


<0.001 


O970 


0.242 


<0.001 


Madrid 


-0.941 


0.275 


0.001 


-1.010 


0.309 


0.001 


Logrono 


-0.333 


0.287 


0.246 


-0.405 


0.287 


0.157 


Majorca 


0.741 


0.246 


0.003 


0.672 


0.263 


O011 


Las Palmas 


-0.333 


0.345 


0.333 


-0.345 


0.359 


0.336 


Sex 


Female 


Male 


-1.735 


0.962 


0.071 


-1.560 


0.886 


0.078 


Age (years) 


0.010 


0.013 


0.459 


0.014 


0.016 


0373 


Financial problems 


Living comfortably 


Doing alright 


0.498 


0.291 


0.087 


0.454 


0.340 


0.182 


Finding it difficult 


0.545 


0.374 


0.144 


0507 


0.376 


0178 


Finding it very difficult 


0.299 


0.525 


0.569 


0226 


0.652 


0729 


Taking psychotropic drugs 


No 


Yes 


0.609 


0.241 


0.011 


0.544 


0.176 


0.002 


Physical health (SF-12) 


Each point on subscale score 


-0.033 


0.008 


<aooi 


-0.030 


0.007 


<0.001 


IVIental health (SF-12) 


Each point on subscale score 


-0.049 


0.009 


<aooi 


-0.048 


0.007 


<0.001 


Dissatisfaction with paid work 


Satisfied 


Dissatisfied 


2.245 


1.041 


0.031 


2.334 


0.953 


O014 


Very dissatisfied 


2.463 


0.928 


0.008 


2.368 


1.094 


0.030 


Jobless 


2.198 


0.527 


< 0.001 


2.072 


0.783 


0.008 


Dissatisfaction with unpaid work 


Satisfied 


Dissatisfied 


0.064 


0.229 


0.779 


0.086 


0.188 


0.646 


Very dissatisfied 


1.129 


0.345 


0.001 


1.054 


0.275 


<0.001 


Not response 


0.144 


0.274 


0.598 


0186 


0.284 


0511 


Sex*age 


0.027 


0.014 


0.059 


0.025 


0.012 


0.039 


Sex'financlal problems 


Living comfortably 


Doing alright 


-0.157 


0.784 


0.841 


-0.276 


0.724 


0.703 


Finding it difficult 


0.133 


0.98 


0.892 


0.043 


0.793 


0.956 


Finding it very difficult 


4.049 


2.231 


0.070 


3.754 


1.607 


0.020 


Age*paid work 


Satisfied 


Dissatisfied 


-0.056 


0.026 


0.035 


-0.059 


0.023 


0.009 


Very dissatisfied 


-0.037 


0.023 


0.102 


-0.035 


0.025 


0.171 


Unemployed 


-0.043 


0.013 


0.001 


-0.042 


0.017 


0.014 


S.E.: Standard error. (*) Multi-level logistic regression with health center as a 


andom component. 











doi:l 0.1 371 /journal.pone.01 06370.t006 



The interaction age*dissatisfaction at work also had an antago- efficient coping strategies for work problems. The fact that anxiety 
nistic effect (Figure 4), a possible reason for which is that younger and depression share many risk factors could be explained in part 
persons have less work experience and may therefore have less by their frequent co-occurrence, although it is also possible that 
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Table 7. Examples of predicted probabilities of anxiety at 
baseline. 



Risk 

score Risk profile 

2.9% A 61 -year-old man living in Madrid 

Very dissatisfied with paid work and satisfied with unpaid work 
*SF-12 mental scale score 61.1 
SF-12 physical scale score 32.3 
Doing alright (financial problems) 

Does not take medication for anxiety, depression or stress 
12.9% A 22-year-old man living in Granada 

Jobless and dissatisfied with unpaid work 
SF-12 mental scale score 32.9 
SF-12 physical scale score 45.8 
Doing alright (financial problems) 

Does not take medication for anxiety, depression or stress 
34.6% A 40-year-old woman living in Las Palmas 

Jobless and satisfied with unpaid work 

SF-12 mental scale score 38.2 

SF-12 physical scale score 26.8 

Doing alright (financial problems) 

Taking medication for anxiety, depression or stress 
84.1% A 37-year-old woman living in Saragossa 

Very dissatisfied with paid work and dissatisfied with unpaid work 

SF-12 mental scale score 23.8 

SF-12 physical scale score 40.2 

Doing alright (financial problems) 

Taking medication for anxiety, depression or stress 

(*)The possible range of the SF-12 is 0-100. High scores indicate good health/ 
well-being. Mean (Standard Deviation) Short Form 12 (SF-12) mental and 
physical subscale scores for Spain are 47.1 (12.4) and 43.8 (11.4), respectively 
[39]. 

doi:l 0.1 371 /journal.pone.01 a6370.ta07 

both are expressions of a latent pathological process [59] . Their 
different phenotypes may be associated with the interaction 
between different risk factors in each case, including those not 
measured in this study (e.g. genes, personality, coping). 

The predicted probability of the onset of anxiety syndromes at 
12 months in an individual could be calculated from the equation 
in Table 3, adding the constant and the corresponding shrunk 
coefficients for each patient and then inverting the "logit" of the 
sum obtained. As these calculations are relatively complex, the 
most reasonable option is a spreadsheet or a web-based calculator. 
Our web-based calculator is available for the onset of major 
depression at 12 months [39], and on the same website (http:// 
www.rediapp.org/predict/Index.php), a calculator will be avaH- 
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able in the coming months for the onset of anxiety syndromes at 
12 months. 

Our results do not address how the predictA-Spain algorithm 
might be implemented in primary care, though this aspect can be 
studied in future research. Recognition by FPs of those patients 
with a greater overall risk (quantitative risk information) as well as 
their risk factors (qualitative risk information) may lead, as in 
cardiovascular disease, to the development of interventions 
tailored for intensity (level) and specificity (profile) of risk. The 
FPs could inform patients about their risk and provide tailored 
counseling and support, thus increasing the patients' empower- 
ment and self-efficacy perception to prevent anxiety syndromes 
[10]. 

The choice of a predicted probability cut-point for making 
clinical decisions depends mosdy on three factors: 1) the validity of 
each cut-point; 2) the available evidence on the effectiveness of 
interventions to prevent anxiety and whether this effectiveness is 
greater for a given predicted probability cut-point; and 3) the 
available evidence on the consequences of false positives and 
negatives regarding patients' health and quality of life and costs for 
patients, health services and society. For example, if an effective 
intervention to prevent anxiety is implemented by FPs, since FPs 
are usually very busy a more specffic cut-point to intervene might 
be preferable (e.g., a predicted probability of 13-14% or greater); 
although this decision should be reconsidered in light of the 
available data on the consequences to health and cost of false 
negative predictions. Nonetheless, if we provide a low cost 
intervention to prevent depression (e.g., an internet-based guided 
self-help) and the consequences linked to the false positives are 
acceptable, a more sensitive cut-point might be interesting (e.g., a 
predicted probability of 10%). Interventions to prevent anxiety 
disorders have been developed mainly in children and adolescents 
[60,61], and less frequently in the elderly [62,63]. However, 
further evaluations are needed of the effectiveness of interventions 
to prevent anxiety disorders in general adult populations. 

Conclusion 

The predictA-Spain risk algorithm is valid to predict anxiety 
syndromes at 12 months. Although external validation is required, 
the predictA-Spain is available for use as a predictive tool in the 
prevention of anxiety syndromes in primary care. 

Acknowledgments 

The authors thank the Primary Care District of Malaga, particularly Jose 
Miguel Morales, Javier Navarro and Maximiliano Vilaseca for their 
support. Finally, we also thank all the interviewers, physicians and patients 
for their participation in the predictA-Spain Study. 

Author Contributions 

Conceived and designed the experiments; JAB FTG CMF EM MJGdGB 
MSG MADE GV JDL PMP MK IN LM GMB. Performed the 
experiments: JAB FTG GMF EM MJGdGB MSG MADB GV PMP. 
Analyzed the data: JAB PMP LM JDL MK IN. Gontributed reagents/ 
materials/analysis tools: None. Wrote the paper: JAB PMP. 



3. Serrano-Bianco A, Palao D, Luciano JV, Pinto-Meza A, Lujan L, ct al. (2010) 
Prevalence of mental disorders in primary care: results from the diagnosis and 
treatment of mental disorders in primary care study (DASMAP). Soc Psychiatry 
Psychiatr Epidemiol 45: 201-210. 

4. Mcrgl R, Seidschek I, Allgaier AK, MoUcr HJ, Hegcrl U, ct al. (2007) 
Depressive, anxiety and somatoform disorders in primary care: prevalence and 
recognition. Depress Anxiety 24: 185—195. 



PLOS ONE I www.plosone.org 



12 



September 2014 | Volume 9 | Issue 9 | el 06370 



Predicting the Onset of Anxiety Syndromes 



5. Fernandez A, BcUon JA, Pinlo-Mcza A, Luciano J V, Autoncll J, ct al. (2010) 
Burden of chronic physical conditions and mental disorders in primary care. 
BrJ Psychiatry 196: 302-309. 

6. Barlow DH, Lehman CL (1996) Advances in the psychosocial treatment of 
anxiety disorders: ImpKcations for national health care. Arch Gen Psychiatry 53: 
727-735. 

7. Kessler RC, Zhao S, Katz SJ, Kouzis AC, Frank RG, et al. (1999) Past-year use 
of outpatient scr\'iccs for psychiatric problems in the national comorbidity 
survey. Am J Psvehiatrv 156: 115-123. 

8. Fernandez A, Haro JM, Martinez-Alonso M, Demyttenaere K, Brugha TS, 
et al. (2007) Treatment adequacy for anxiety and depressive disorders in six 
European countries. BrJ Psychiatry 190: 172—173. 

9. Schmidt NB, Zvolensky MJ (2007) Risk factor research and prevention for 
anxiety disorders: introduction to the special series on risk and prevention of 
anxiety pathology. Behav Modif 31: 3-7. 

10. BellonJA, Conejo-Ceron S, Moreno-Peral P, King M, Nazareth I, et al. (2013) 
Preventing the onset of major depression based on the level and profile of risk of 
primary care attendees: protocol of a cluster randomised trial (the prcdictD- 
CCRT study). BMC Psychiatry 13: 171. 

11. King M, Bottomley C, BellonJA, Torres-Gonzalez F, Svab I, et al. (2011) An 
international risk prediction algorithm for the onset of generalized anxiety and 
panic syndromes in general practice attendees: prcdictA. Psychol Med 41: 1625— 
1639. 

12. Bellon JA, Moreno-Kiistner B, Torres-Gonzalez F, Monton-Franco C, 
GildeGomez-Baxrag^ MJ, et al. (2008) Predicting the onset and persistence 
of episodes of depression in primary care: the PREDICT-Spain Study. 
Methodology. BMC Pubhc Health 8: 256. 

13. Spitzer RL, Kroenke K, Williams JBW (1999) Validation and Utility of a Self- 
report Version of PRIME-MD. The PHQ_ Primary Care Study. JAMA 282: 
1737-1744. 

14. Baca E, Saiz J, Agiiera L, Caballero L, Fernandez-Liria A, et al. (1999) 
Validacion de la version espaiiola del PRlME-MD: un procedimiento para el 
diagnostico de trastornos mentales en atencion primaria. Actas Esp Psiquiatr 27: 
375-383. 

15. Karasek RA, Theorell T (1990) Healthy Work: Stress, Productivity, and the 
Reconstruction of Working Life. New York, NY: Basic Books. 

16. Welch S, Lewis G (1998) Poverty, unemployment, and common mental 
disorders: population based cohort study. BMJ 317: 115-119. 

17. Jenkinson C, Layte R, Jenkinson D, Lawrence K, Petersen S, et al. (1997) A 
shorter form ealth survey: can the SF-12 replicate results from the SF-36 in 
longitudinal studies? J Public Health Med 19: 179-186. 

18. Gandek B, Ware JE, Aaronson NK, Apolone Ci, BjornerJB, et al. (1998) Cross- 
validation of item selection and scoring for the SF-12. Health Sur\'ey in nine 
countries: results from the IQpLA Project. International Quality of Life 
Assessment.J Clin Epidemiol 51: 1171-1178. 

19. Barbor TF, de la Fuente JR, Saunders J, Grant M (1989) The Alcohol Use 
Disorders Identification Test: Guidelines for the Use in Primary Health Care. 
Genova: World Health Organization. 

20. Rubio Valladolid G, Bermejo Vicedo J, Caballero Sanchez-Serrano MC, Santo- 
Domingo CarrascoJ (1998) Validation of the alcohol use disorders identification 
test (AUDIT) in primary care. Rev Clin Esp 198: 1 1-14. 

21. Perula-de 'Forres LA, Fernandez-CiareiaJA, Arias- Vega R, Muriel- Palomino M, 
Marquez-ReboUo E, et al. (2005) Validity of AUDIT test lor detection of 
disorders related with alcohol consumption in women. Med Clin (Bare) 125: 
727-730. 

22. Arroll B, Khin N, Kerse N (2003) Screening for depression in primairy care with 
two verbally asked questions: cross sectional study. BMJ 327: 1144-1146. 

23. WHO. Composite International Diagnostic Interview (CIDI) (1997) Version 2.1. 
Geneva: World Health Organization. 

24. Reynolds CF, Frank E, Thase ME, Houck PR, Jennings JR, et al. (1988) 
Assessment of sexual function in depressed, impotent, and healthy men: factor 
analysis of a Brief Sexual Function Questionnaire for men. Psychiatry Res 24: 
231-250. 

25. Robins LN, Wing J, Wittchen HU, Helzcr JE, Babor TF, et al. (1988) The 
Composite International Diagnostic Interview. An epidemiologic instrument 
suitable for use in conjunction with different diagnostic systems and in different 
cultures. Arch Gen Psychiatry 45: 1069-1077. 

26. Rubio-Stipec M, Bravo M, Canino G (1991) The Composite International 
Diagnostic Inter\dew (CIDI): an epidemiologic instrument suitable for using in 
conjunction with different diagnostic systems in different cultures. Acta Psiquiatr 
Psieol Am Eat 37: 191-204. 

27. Fink LA, Bernstein D, Handelsman L, Foote J, Lovejoy M (1995) Initial 
reliability and validity of the childhood trauma interview: a new multidimen- 
sional measure of childhood interpersonal trauma. Am J Psychiatry 152: 1329— 
1335. 

28. King M, Speck P, Thomas A (1995) The Royal Free interview for religious and 
spiritual beliefs: development and standardization. Psychol Med 25: 1 125-1 134. 

29. Tyrer P (1990) Personality disorder and socIeJ functioning. In Peck DF, Shapiro 
CM, editors. Measuring Human Problems: A Practical Guide. Chichester: 
Wiley & Sons. 119-142 p. 

30. Qureshi N, Bethea J, Modell B, Brerman P, Papageorgiou A, et al. (2005) 
Collecting genetic information in primary care: evaluating a new farrnly history 
tool. Fam Pract 22: 663-669. 



31. Sproston K, Primatesta P (2003) Health Sur\-ey lor England 2002: a Survey 
Carried out on Behalf of the Department of Health. Volume 1: The Health of 
Children and Young People. London, England: The Stationery Office. 

32. Brugha T, Bebbington P, Termant C, Hurry J (1985) The List of Threatening 
Experiences: a subset of 12 Ufe event categories with considerable long-term 
contextual threat. Psychol Med 15: 189-194. 

33. Motrico E, Moreno-Kustner B, Luna JD, Torres-Gonzalez F, King M, et al. 
(20 1 3) Psychometric properties of the List of Threatening Experiences-LTE and 
its association with psychosocial factors and mental disorders according to the 
different scoring methods. J Affect Disord 150: 931-940. 

34. Janssen 1, Hanssen M, Bak M, BijI RV, de Graaf R, et al. (2003) Discrimination 
and delusional ideation. BrJ Psychiatry 182: 71-76. 

35. Blaxter M (1990) Health and Lifestyles. London, England: Roudedge. 

36. Heatherton TF, Kozlowski LT, Freeker RG, Riekert W, Robinson J (1989) 
Measuring the heaviness of smoking: Using self-reported time to the first 
cigarette of the day and number of cigarettes per day. BrJ Addict 84: 791-799. 

37. Greenland S (1989) Modeling variables selection in epidemiologic analysis. 
AmJ Public Health 79: 340-349. 

38. Bebbington P, Dunn G, Jenkins R, Lewis G, Brugha T, et al. (2003) The 
influence of age and sex on the prevalence of depressive conditions: report from 
the National Survey of Psychiatric Morbidity. Int Rev Psychiatry 15: 74-83. 

39. BellonJA, Luna J, King M, Moreno-Kustner B, Nazareth I, et al. (2011) 
Predicting the onset of major depression in primary care: international validation 
of a risk prediction algorithm from Spain. Psychol Med 41: 2075-2088. 

40. Hinz A, KittelJ, Karoff M, Sehwarz R (2004) Age and sex dependencies of 
anxiety and depression in cardiologic patients compared with the general 
population. Psychosoc Med 1: Doc09. 

41. Heman MA, Hernandez-Diaz S, Robins JM (2004) A structural approach to 
selection bias. Epidemiology 15: 615-625. 

42. BeU6n JA, Luna JD, Moreno-Kiistner B, Monton-Franco C, GildeGomez- 
Barragan MJ, et al. (2010) Psychosocial and sociodemographie predictors of 
attrition in a longitudinal study of major depression in primary care: The 
predictD-Spain study. J Epidemiol Community Health 64: 874-884. 

43. Harrell EE (2001) Regression Modelling Strategies. New York, NY: Springer. 

44. Copas JB (1983) Regression, prediction and shrinkage. J R Statist Soe B 45: 
311-354. 

45. Cooper H, Hedges LV (1994) The Handbook of Research Synthesis. New York, 
NY: Russell Sage Foundation. 

46. StataCorp (2011) Stata statistical software: Release 12. In Station C, editor. 
Texas: StataCoip LP. 

47. Snijders TAB, Bosker RJ (1999) Multilevel Analysis: An Introduction to Basic 
and Advanced Multilevel Modelling. London: Sage Publications. 

48. Altman DG, Royston P (2000) What do we mean by validating a prognostic 
model? Stat Med 19: 453-473. 

49. Lee ML, Yano EM, Wang M, Simon BF, Rubenstein LV (2002) What patient 
population does visit-based sampling in primary care setting represent? Med 
Care 40: 761-770. 

50. Wyke S, Hunt K, Walker J, Wilson P (2003) Frequent attendance, 
sociocronomie status and burden of ill health. An investigation in the west of 
Scotland. Eur J Gen Pract 9: 48-55. 

5 1 . King M, Walker C, Levy G, Bottomley C, Royston P, et al. (2008) Development 
and vahdation of an international risk prediction algorithm for episodes of major 
depression in general practice attendees: the PredictD study. Arch Gen 
Psychiatry 65: 1368-1376. 

52. Conroy RM, Pyorala K, Fitzgerald AP, Sans S, Menotti A, et al. (2003) 
Estimation of ten-year risk of fatal cardiovascular disease in Europe: the SCORE 
project. Eur Heart J 24: 987 1003. 

53. Moons KG, Royston P, \7Tgou\ve Y, Grobbee DE, Altman DG (2009) 
Prognosis and prognostic research: what, why, and how? BMJ 338: b375. 

54. Moreno-Peral P, Conejo-Ceron S, Motrico E, Rodriguez-Morejon A, 
Fernandez A, et al. (2014) Risk factors for the onset of panic and generalized 
anxiety disorders in the general adult population. A systematic review of cohort 
sttidies. J Affect Disord 168: 337-348. 

55. Alonso J, Angermeyer MC, Bemert S, Bruffaerts R, Brugha TS, et al. (2004) 
Psychotropic drug utilization in Europe: results from the European Study of the 
Epidemiology of Mental Disorders (ESEMeD) project. Acta Psychiatr Scand 
SuppI 420: 55-64. 

56. European Commission EUROSTAT (2013) Income distribution statistics, 
http:/ / epp.eurostat.ee. europa.eu/ statisties_explained/index.php/Income_ 
distribution_statisties. Accessed 1 October 2013. 

57. Keyes CLM, Dhingra SS, Simoes EJ (2010) Change in Level of Positive Mental 
Health as a Predictor of Future Risk of Mental Illness. AmJ Public Health 100: 
2366-2371. 

58. King M, Nazaretii I, Levy G, Walker C, Morris R, et al. (2008) Prevalence of 
common mental disorders in general practice attendees across Europe. 
BrJ Psychiatry 192: 362-367. 

59. Mennin DS, Heimberg RG, Fresco DM, Rittcr MR (2008) Is generalized 
anxiety disorder an anxiety or mood disorder? Considering multiple factors as 
we ponder the fate of GAD. Depress Anxiety 25: 289-299. 

60. Christensen H, Pallister E, Smale S, Hickie IB, Calear AL (2010) Community- 
based prevention programs for anxiety and depression in youth: a systematic 
review. J Prim Prev 31: 139-170. 



PLOS ONE I www.plosone.org 



13 



September 2014 | Volume 9 | Issue 9 | e106370 



Predicting the Onset of Anxiety Syndromes 



61. Calear AL, Christensen MJA (2010) Review of intcrnct-bascd prevention and 
treatment programs for anxiety and depression in ehildrcn and adolescents. 
MedJ Aust 192: S12-S14. 

62. Van't Veer-Tazelaar PJ, van Marwijk HWJ, van CJppen P, van Hout HP, van 
der Horst HE, et al. (2009) Stepped-Care prevention of anxiety and depression 
in late life. Arch Gen Psychiatry 66: 297-304. 



63. Van't Veer-Tazelaar PJ, van Alarwijk H\\(J, van Oppen P, van der Horst HE, 
Smit F, ct al. (2011) Anxiety and Depression Has Sustained Elfects Over 24 
Months: A Pragmatic Randomized Trial. Am J Geriatr Psychiatry 19: 230-239. 



PLOS ONE I www.plosone.org 



14 



September 2014 | Volume 9 | Issue 9 | el 06370 



